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ABSTRACT 

We present chemical abundance measurements from high resolution observations of 7 sub- 
damped Lyman-a absorbers and 1 damped Lyman-a system at z < 1.5. Three of these objects 
have high metallicity, with near or super-solar Zn abundance. Grids of Cloudy models for 
each system were constructed to look for possible ionization effects in these systems. For 
the systems in which we could constrain the ionization parameter, we find that the ionization 
corrections as predicted by the Cloudy models are generally small and within the typical error 
bars (~ 0.15 dex), in general agreement with previous studies. The Al III to Al II ratio for 
these and other absorbers from the literature are compared, and we find that while the sub- 
DLAs have a larger scatter in the Al III to Al II ratios than the DLAs, there appears to be 
little correlation between the ratio and Nh i- The relationship between the metallicity and the 
velocity width of the profile for these systems is investigated. We show that the sub-DLAs that 
have been observed to date follow a similar trend as DLA absorbers, with the more metal rich 
systems exhibiting large velocity widths. We also find that the systems at the upper edge of this 
relationship with high metallicities and large velocity widths are more likely to be sub-DLAs 
than DLA absorbers, perhaps implying that the sub-DLA absorbers are more representative 
of massive galaxies. 
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1 INTRODUCTION 

The chemical composition of galaxies provides important clues 
into galaxy formation and evolution. Absorption systems in QSO 
spectra provide an opportunity to study the interstellar medium 
(ISM) and therefore the chemical enrichment of galaxies. Damped 
Lyman-a (DLA) and sub-Damped Lyman-a (sub-DLA) absorbers 
are 2 classes of absorbers seen in the spectra of background QSOs. 
These high neutral hydrogen column density absorbers, particularly 
DLA s with log Nh i > 2 0.3 and sub-DLAs with 19 < log N H i < 
20.3 jPeroux et alj|200ll) are expected to contain the majority of 
the neutral gas in the universe, while the majority of the baryons 
are thought to lie in the highl y ionized and d iffuse Lyman-a forest 
clouds with log N H i < 14 JPetitiean et al]|l993l) . The DLA and 
sub-DLA systems are believed to be strongly related to galaxies 
at any redshift. Indeed, several DLA host galaxies have been con- 



firmed through deep imagin g at z < 0.5 eg; dChen and LanzettU 
l2003l : lGharanfoli et al.ll2007t) . 

DLA systems have long been the preferred class of 
absorbers for chemical abunda nce studies. Observations by 
Dessauges-Zavadskv et al. I d2003h of 12 sub-DLA absorbers at 
1.8 < z < 4.3 showed that these systems also contain metal 
line absorption featu res. Based on Fe II abundance measurements, 
IPeroux et al. (2003a) suggested stronger metallicity evolution in 
sub-DLA systems than DLA systems. To date, very few observa- 
tions have been made of z < 1.5 sub-DLAs. Redshifts z < 1.5 
span 70% of the age of the universe (assuming a concordance 
cosmology of fl m = 0.3, Qa = 0.7). Clearly this redshift range 
is important to understanding the nature of sub-DLA systems and 
galactic chemical evolution as well. 

Ionization effects can alter the abundances inferred for QSO 
absorbers. For DLAs, it is typically assumed that the ionization 
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corrections are small due to their large H I column densities that 
shield the cloud from ionizing photons. Dessauges-Zavadsk v et al.l 
d2003h showed that the ionization corrections for sub-DLA systems 
are also typically s mall, with < 0.2 dex corrections to most ions. It 
was also shown bv lPerouxetai] (2003b) that the contributions from 
DLAs to the neutral gas budget of the universe at z > 3.5 was lower 
than at z < 3.5 and that sub - DLAs harbored thi s remaining HI. 



Recently 
(2007); P 



Recently Peroux et al. (2006a); York et al. (2006); Kulkarni et al. 



; Prochaska et al. ( 2006); Khare et al. (2006) have suggested 
that sub-DLAs may be more metal rich than DLAs, and may con- 
tribute significantly to the cosmic metal budget. Sub-DLAs have 
thus become a topic of much interest. 

Many elements have been detected in QSO absorber systems, 
including C, N, O, Mg, Si, S, Ca, Ti, Cr, Mn, Fe, Ni, and Zn. Zn is 
often the preferred tracer of the gas-phase metallicity as it is rela- 
tively undepleted in the Galactic ISM, especially when the fraction 
of H in molecular form is low which is true in most DLAs. Zn also 
tracks the Fe abundance in Galactic stars, and the lines of Zn II 
AA 2026,2062 are relatively weak and typically unsaturated. These 
lines can also be covered with ground based spectroscopy over a 
wide range of redshifts, from 0.65 < z < 3.5, which covers a large 
portion of the history of the universe. Relative abundances of re- 
fractory elements such as Cr and Fe relative to Zn also relate to the 
amount of extinction. 

As galaxies mature, the ongoing processes associated with 
star formation such as stellar winds, supernovae explosions etc. ex- 
pel metals into the interstellar medium. Over time, the interstellar 
metallicity is expected to reach a near solar value. Several models 
based on both analytical and semi-analytical calculations, as well as 
hydrodynamical simulations predict this increase to approximately 
solar values. Recent findings indicate that the Nh i weighted mean 
metallicity (Z) of DLAs is well below solar (~ 0. IZq at z — 2.0 
and ~0.16Zq at z = 1) and shows little if any evolution that is 
predicted by models dKulkarniet ai1 l2005). Some lines of investi- 
gation are now suggesting that sub-DLAs may contain a significant 
portion of these "missing metals", and may be more representa- 
tive of typical m a ssive spirals o r ellipticals I Kulkarn i et al.ll2007l ; 
iKhare et al.l2006h . lWolfefdl993l) claimed that the DLAs can be con- 
sidered as progenitors to normal modern spirals, and argued that the 
line profiles suggest rotation curves (in absorption). Others how- 
ever note that the absence of large galaxies in a number of searches 
for DLA host galaxies, suggesting that many are dwarf galaxies 
often with low impact parameters and lost in the l i ght of the QSO 



men witn low impact parameters ana lost m the lignt ot tne QaU 
York et alJl986l ; lHopkins. Rao, fc Turnshekl2005hlKulkarni et all 



20071) how that based on available spectroscopic measurements, 
at z < 1 the contribution of sub-DLAs to the total metal budget 
may be several times greater than that of DLAs. For a constant 
dust-to-gas ratio, the dust obscuration bias is also expected to be 
less for sub-DLAs than for DLAs, therefore one might expect to be 
able to find metal rich sub-DLAs more frequently if dust obscura- 
tion effects are important. 

Although difficult to constrain, it is generally acknowledged 
that there may be a significant ionization fraction in the sub- 
DLA systems, especially for systems with lower Nh i dPeroux et al.l 
l2003al ; IO'Meara et all 2006). It has also been claimed that the sub- 
DLAs contain ~15% of the H atoms in the universe, and that they 
may be a d istinct population of absorbers separated from the DLA 
absorbers dO'Meara et alj|2006t) . The H I column density distribu- 
tion for QSO absorbers shows that systems with 19 < log Nh i < 
20.3 are ~8 times more numerous than the higher column den- 
sity classical damped Lya systems, making sub-DLAS more read- 
ily available probes of neutral gas in the distant Universe. 



In this paper we present high-resolution measurements of 7 
sub-DLAs and 1 DLA taken with the Magellan Inamori Kyocera 
Echelle (MIKE) spectrograph on the 6.5m Clay telescope at Las 
Campanas observatory. The structure of this paper is as follows. 
In § 2, we discuss details of our observations and data reduction 
techniques. In § 3, column density determinations are discussed. § 4 
gives information on the individual objects, while in § 5 we discuss 
the abundances of these absorbers, the Al III to Al II ratio, and 
ionization corrections, in § 6 we draw some general conclusions. 



2 OBSERVATIONS AND DATA REDUCTION 

The spectra presented here were obtained over 2 separate epochs, 
2005 Sep. and 2006 February, respectively, with the 6.5m Mag- 
ellan Clay telescope at Las Campanas Observatory. These objects 
were observed wi th the Magellan Inam ori Kyocera Echelle spectro- 
graph (MIKE) dBernstein et alj|20030 . This is a double sided spec- 
trograph with both a blue and a red camera, providing for simulta- 
neous wavelength coverage from ~3340 A to ~9400 A . Targets 
were observed in multiple exposures of 1800 to 2700 sec to mini- 
mize cosmic ray defects. The seeing was typically < 1", averaging 
~ 0.7". All of the target QSOs were observed with the l"x5" slit 
and the spectra were binned 2x3 (spatial by spectral) during read- 
out. The resolving power of the MIKE spectrograph is ~ 19,000 
and ~25,000 on the red and blue sides respectively with a l"x5" 
slit. Table 1 gives a summary of the observations. 

These spectra were reduced using the MIKE pipeline reduc- 
tion code in IDL developed by S. Buries, J. X. Prochaska, and R. 
Bernstein. Wavelengths were calibrated using a Th-Ar comparison 
lamp taken after each exposure. The MIKE software first bias sub- 
tracts from the overscan region and flat-fields the data, then sky- 
subtracts and extracts the spectral orders using the traces from flat 
field images. The pipeline calibration software also corrects for he- 
liocentric velocities and converts the wavelengths to vacuum val- 
ues. Each individual order was then extracted from the IDL struc- 
ture created by the pipeline software and combined in IRAF using 
rejection parameters to reduce the effects of cosmic rays. These 
combined spectra were then normalized using a Legendre polyno- 
mial function to fit the continuum. Typically, these functions were 
of order five or less. 

Our sample consists of 7 sub-DLAs and 1 DLA at z<1.5, 
where there are few existing abundance measurements, including 3 
absorbers with z<1.0 where there are even fewer a bu ndance mea- 
surem ents dKhare et ai]|2004l; IPeroux et aT1l2006dlbl iPettini et all 
2000). The objects from our sample were chosen partly because 
of strong lines of Mg II AA 2796, 2803, or Fe II A 2600 seen 
in SPSS spectra, or otherwise previously measured and given in 
iRao. Turnshek. fc Nestoil j20 06). These sub-DLAs and the two ob- 
served in IPeroux et alj d2006al lbh comprise all of high resolution 
observations of sub-DLAs at z < 1.5 to date. Four of the tar- 
get QSOs in this sample were observed in the Sloan Digital Sky 
Survey (SDSS) and all have a kn own Nh i from HST spectra. H I 
column densities were taken from Rao, Turnshe k. fc Nestoil (2006) 
and references therein. We provide an appendix at the end of this 
paper showing the fits of the Lyman-Q profiles from the UV spectra. 
Throughout this paper, the QSO names are given in 12000 coordi- 
nates. 
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Table 1. Summary of Observations 



QSO 
J2000 


Original or SDSS ID 


RA 


Dec 


my 


z em 




Nhi 

cm" 2 


Exposure Time 
Sec 


Hpoch 


Mg 11 Reference 


Q0354-2724 


Q0352-272 


03:54:05.9 


-27:24:25.7 


17.9 


2.823 


1.4051 


20.18±0.15 


10800 


2005 Sep, 2006 Feb. 


1 


Q0826-2230 


Q0823-223 


08:26:01.5 


-22:30:26.2 


16.2 


>0.911 


0.9110 


19.04±0.04 


8627 


2006 Feb. 


2 


Q 1009-0026 


SDSS 1100930.46-002619.1 


10:09:30.4 


-00:26:19.1 


17.4 


1.244 


0.8426 


20.20±0.06 


8100 


2006 Feb. 


3 














0.8866 


19.48±0.05 






3 


Q1010-0047 


SDSS J101033. 44-004724.5 


10:10:33.4 


-00:47:24.5 


18.0 


1.671 


1.3270 


19.81±0.05 


7191 


2006 Feb. 


3 


Q 1224+0037 


SDSS J122414.29+003709.0 


12:24:14.3 


00:37:09.0 


18.7 


1.482 


1.2346 


20.88±0.05 


5400 


2006 Feb. 


3 














1.2665 


20.00±0.07 






3 


Q2331+0038 


SDSS J233121. 81+003807 .4 


23:31:21.8 


00:38:07.4 


I7.fi 


1.486 


1.1414 


20.00±0.05 


80(10 


2005 Sep 


3 



Notes - Emission redsh ifts are from the reference below , or if t he t arget was observ ed i n the SDSS, the emission redshift is from lSchneider et alj feOOSl) 
Mg II References. - m lSargent, Steidel, fc Boksenberd fl989l) . (2) lFalomolfl990l) . ffl lNestorlfeoWl) 



Table 2. Rest-frame equivalent width measurements of key metal lines. Values and lcr errors are in units of mA . 



QSO 


z af>s 


Mgl 


Mg II 


Mgll 


Al 11 


Al III 


Aim 


Si 11 


Si 11 


Call 


Call 


Cr 11 






2852 


2796 


2803 


1670 


1854 


1862 


1526 


1808 


3933 


3969 


2056 


Q0354-2724 


1.4051 


537 ± 9 


2665 ± 14 


2394 ± 10 


< 1114" 


< 412 a 


< 607 a 


< 1670" 








64 ± 8 


Q0826-2230 


0.9110 


123 ± 9 


1280 ± 16 


892 ± 15 




<9 


<9 




<10 


22 ± 4 


14 ± 4 


<5 


Q 1009-0026 


0.8426 


77 ± 14 


713 ± 13 


551 ± 13 




77 ± 8 


33 ± 6 










<10 


01009-0026 


0.8866 


308 ± 22 


1792 ± 16 


1508 ± 18 




146 ± 10 


88 ± 13 




<12 


103 ± 16 


83 ± 14 


<5 


Q1010-0047 


1.3270 


343 ± 18 


2040 ± 14 


1768 ± 15 


883 ± 17 


279 ± 9 


141 ± 8 


681 ± 21 


44 ±11 






<9 


Q 1224+0037 


1.2346 


212 ± 15 


1086 ± 15 


1032 ± 18 


513 ± 17 


159 ± 10 


99 ± 9 


477 ± 22 


52 ± 7 






47 ± 7 


Q1224+0037 


1.2665 


188 ± 47 


2181 ± 28 


1876 ± 25 


827 ± 15 


84 ±24 


<10 


<1140 6 


< 12 






<9 


Q2331+0038 


1.1414 


299 ± 19 


2534 ± 43 


2066 ± 47 


424 ± 85 


105 ± 50 


52 ± 25 




<13 


257 ± 23 


141 ± 26 


<9 


QSO 


z af,s 


Mn 11 


Mnll 


Mnll 


Fen 


Fell 


Fell 


Fell 


Fc 11 


Fe [1 


ZnH c 


Znn d 






2576 


2594 


2606 


2260 


2344 


2374 


2382 


2586 


2600 


2026 


2062 


Q0354-2724 


1.4051 


126 ± 13 


100 ± 10 


75 ± 8 


109 ± 9 


1249 ± 18 


739 ± 17 


1657 ± 21 


1 143 ± 16 


1715 ± 18 


84 ± 10 


68 ± 15 


Q0826-2230 


0.9110 


<5 


<5 


<5 


<4 


161 ± 21 


39 ± 15 


444 ± 14 


102 ± 10 


353 ± 20 


44 ± 8 


26 ± 14 


Q 1009-0026 


0.8426 


27 ± 4 


42 ± 6 


<8 


31 ± 6 


275 ± 6 


179 ± 5 


389 ±4 


265 ± 6 


403 ± 5 


<13 


<10 


Q 1009-0026 


0.8866 


<5 


<5 


<5 


23 ±4 


625 ± 10 


255 ± 10 


1000 ± 12 


568 ± 1 1 


981 ± 18 


41 ± 6 


<5 


Q1010-0047 


1.3270 


<15 


<15 


<15 





905 ± 32 


396 ± 24 


1230 ± 14 


741 ± 17 


1269 ± 13 


<9 


<6 


Q 1224+0037 


1.2346 


<29 


<29 


<29 


<40 


632 ± 34 


429 ± 35 


758 ± 24 


658 ± 29 


827 ± 29 


<14 


30 ± 7 


Q1224+0037 


1.2665 


<31 


<31 


<31 


<30 


757 ± 28 


271 ±40 


1206 ± 36 


736 ± 65 


1385 ± 27 


<13 


<9 


Q233 1+0038 


1.1414 


<13 


<13 


<13 


< 12 


509 ± 13 


324 ± 17 


1206 ± 46 


662 ± 36 


1285 ± 36 


29 ± 10 


<7 



"Lines were blended with Lyman-a forest lines. 
b Blended with another feature. 

c This line is a blend with Mg IA 2026, although the Mg I contribution is judged to be insignificant in all cases. 
d As this line is blended with the Cr II A 2062 line, this value represents the total EW of the line. 



3 DETERMINATION OF COLUMN DENSITIES 

Column densities were determined from profile fitting with the 
package FITS6P dWeltv et al.lll99lh . which has evolved from the 
code by Vidal-Madjar et al (1977). FITS6P iteratively minimizes 
the x 2 value between the data and a theoretical Voigt profile that is 
convolved with the instrumental profile. The profile fit used multi- 
ple components, tailored to the individual system. For the central, 
core components, the Doppler parameters (6 e //) and radial veloci- 
ties were determined from the weak and unsaturated lines, typically 
the Mg I A 2852 line. For the weaker components at higher radial 
velocities the 6 e // and component velocity values were determined 
from stronger transitions such as the Fe II AA 2344, 2382 lines and 
the Mg II AA 2796, 2803 lines. A set of b e // and v values were thus 
determined that reasonably fit all of the lines observed in the sys- 
tem. Th e atomic data u sed in line identification and profile fitting 
are from lMortonl fc003h . 

In general, if a multiplet was observed, the lines were fit si- 
multaneously until convergence. For all of the systems, the Fe II A 



2344, 2374, 2382 lines were fit simultaneously to determine a set of 
column densities that fit the spectra reasonably well. Similarly, the 
Mg II AA 2796, 2803 lines were also fit together. The Zn II A 2026 
line is blended with a line of Mg I. The Mg I contribution to the 
line was estimated using the Mg I A 2852 line, for which fA ~32 
times that of the A Mg 1 2026 line. The Zn II components were then 
allowed to vary while the Mg I components were held fixed. No n 
was determined by simultaneously fitting the Cr II A 2056 line and 
the blended Cr II + Zn II A 2062 line, where the contribution from 
Zn II was estimat ed from the Zn II + Mg I A 2026 line. See also 
iKhare et"al] d2004l) for a discussion of the profile fitting scheme. In 
this paper we adopt the standard notation 

[X/Y] = logiNx/Nn i) - log(X/H) Q . 

Solar system abundances have been adopted from ?. 

Using the package SPECP, also developed by D.E. Welty, we 
determined equivalent widths as well as column densities via th e 
apparent optical depth method (AOD) dSavage fc Sembachll 1996b . 
We present the rest-frame equivalent widths (Wo) of the lines in 
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Table 4. The la errors for the equivalent widths are given also and 
reflect both uncertainties in the continuum level and in the pho- 
ton noise. If a certain line was not detected, the limiting equivalent 
width was determined from the local signal to noise ratio (S/N). As- 
suming a linear curve of growth, a 3a upper limit was also placed 
on the column density. 



4 DISCUSSION OF INDIVIDUAL OBJECTS 
4.1 Q0354-2724 (z em = 2.823) 

This QSO (Mv=17.9) has a known sub-DLA with lo 
N H i=20.18 at z a63 =1.4051 dSargent, Steidel. fe Boksenberdl 19891 ; 



iRao. Turnshek. fe NestoJ l2006t) . We observed this target during 
both the 2005 Sep. and 2006 Feb. epochs for a total exposure time 
of 10800 sec. We detect strong absorption features of Mg I, Mg II, 
Cr II, Mn II, Fe II, and Zn II. The S/N varies from ~20 at ~ 4900 A 
to ~50 at ~5700 A . Several lines of interest such as Al II A 1670, 
Al m A A 1854, 1862 and Si II A 1808 were inside the Lyman-a for- 
est due to the relatively high redshift of the background QSO. Fig- 
ure 1 shows velocity plots of several lines of interest for this system. 
There are no detectable components from the absorption system in 
any species outside the plotted range. Sargent et al. (1989) reported 
the existense of a C IV absorption system at z a t s = 2.1442 based 
on low-resolution (~ 6 A ) observations. Based on these higher 
resolution spectra, we believe that this system was misidentified. 
There is a small feature at ~350 km s _1 of the Zn II A 2026 line 
(X bs= 4878.7 A ), as can be seen in figure 1, with W b s =67 mA 
. However, if this line is assumed to be the C IV A 1550 line, then 
the A 1548 line should be ~4 A blueward of the Zn II 2026 line. 
At the resolution of Sargent et al. (1989), the Zn II line from the 
Zabs = 1.4051 sub-DLA was likely confused for this C IV A 1548 
line. This feature is unidentified, it does not match any possible 
lines from the other known systems. 

This system shows a complex velocity structure with a total 
of at least 8 components necessary to fit the observed profile. The 
total profile spans ~360 km s~ 1 (see § 5.4 for a more detailed dis- 
cussion of how we derived the velocity width). We detect Zn II at 
S/N ~ 20 with log N Zn n=12.73±0.03 for this sub-DLA giving a 
nearly solar abundance [Zn/H]=-0.08. The Mg I contribution to the 
blended Zn II + Mg I A 2026 line was estimated from the Mg I 
A 2852 line. In theory, if the Mg I A 2852 line is saturated, which 
is rare but does occur occasionally for DLAs, the estimated con- 
tribution to the blended Zn II + Mg I A 2026 line from the Mg 
I A 2852 line could be lower than its true value. The EW of the 
Mg I A 2852 line was W resi = 537 mA , which is below the 600 
mA limit where the line begins to becomes significantly saturated 
dYork et all 2006). This system also shows moderate iron depletion 
with [Fe/Zn]=-0.43. Table 3 gives the results of the profile fitting 
analysis for each component. 



4.2 Q0826-2230 (z em > 0.911) 

This is a BL Lac object. Due to the ambiguous emission redshift 
for this object, it is possible that this sub-DLA system is associated 
with the source itself. This object has a weak sub-DLA with log 
N H i= 19.04 at z =0.91 10 dFalomcll99(j|Rao. Turnshek. fc Nestol 
2006). A total of 7 components were used in the theoretical profile. 
The total velocity width of the profile was ~322 kms _1 .Fi gure 
2 shows velocity plots for some of the lines of interest. We detect 




Figure 1. Velocity plots for several lines of interest in the z = 1 .405 1 system 
in the spectrum of Q0354-2724. The solid green line indicates the theoret- 
ical profile fit to the spectrum, and the dashed red line is the continuum 
level. The vertical dotted lines indicate the positions of the components that 
were used in the fit. In the case of the Zn II AA 2026,2062 lines, the dashed 
vertical lines indicate the expected positions of the components of the Mg I 
(A 2026) and Cr II (A 2062) contributions respectively. 



strong features of Mg II, Fe II, and Ca II. Mg I is also unambigu- 
ously detected, but the majority of the profile is contained in the 
component at v — —65 km s , We also detect Zn II A 2026 and 
2062 at S/N ~ 40 in the region, with the majority of the contribution 
from the v — —65 km s _1 component. There was a null detection 
of the Cr II A 2056 line, which is the intrinsically strongest line 
available, so only a limit could be placed on No n. 

This system has super-solar metallicity with [Zn/H]=+0.68 
(see the appendix for a discussion of the H I values). The weaker 
Fe II AA 2249, 2260 lines were not detected in this system, but the 
Fe II A 2374 line was unsaturated (W re3 t=39 mA , spread over 1 A 
) and an accurate column density could be determined. This system 
appears to have a high depletion with [Fe/Zn]=-1.63. The measured 
Ca II column density of log Nc a n= 1 1 .75 likely underestimates the 
true Ca column density because the ionization potential for Ca II is 
1 1.868 eV is less than the ionization potential of H I. Table 4 gives 
the results of the Voigt profile fitting for each component. 

Grids of Cloudy models (Fer land et alj |j998) were con- 
structed for this system which were tailored to fit the Nh i and 
metallicity observed. The Al II A 1670 line for this system was not 
in the wavelength section covered, and only limits could be placed 
on the Al III A A 1854, 1862 lines. We were however able to con- 
strain the ionization parameter U based on the Mg II to Mg I ratio 
observed. Since the Mg II AA lines were saturated, this ratio is only 
a lower limit, and the subsequent ionization corrections may be dif- 
ferent. However, the Mg II lines were not optically thick at this 
resolution, and the true column density may not be much higher 
than the value given. The AOD column density also agrees very 
well with the profile fitting column density. Based on the model 
predictions, the ionization parameter is small and the ionization 
corrections for Zn are negligible (<0.1 dex) for a large range of 
ionization parameters. See §5.3 for a more thorough account of the 
Cloudy modeling and ionization corrections. 
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Table 3. Column Densities for the z =1.4051 absorber in Q0354-2724. Velocities and b e ff values are in units of km 
s -1 , and column densities are in units of cm -2 . 



Vel 


b eff 


Mgl 


Mgll 


Cr II 


Mnll 


Fell 


Znll 


-151 


10.9 




(3.32±0.53)E+12 


(2.07±0.68)E+12 




(3.24±0.90)E+12 




-133 


5.0 




(3.22±0.77)E+12 










-95 


18.7 


(1.04±0.06)E+12 


>2.85E+13 




(7.09±2.09)E+11 


(4.89±0.20)E+13 


(7.61±1.29)E+11 


-64 


6.4 


(1.63±0.38)E+11 


>1.57E+14 






(9.20±1.10)E+12 




-39 


6.3 


(3.67±0.43)E+11 


>1.21E+14 




(6.10±1.53)E+11 


(5.66±0.39)E+13 







17.4 


(6.64±0.63)E+11 


>4.71E+13 


(2.73±0.79)E+12 




(4.62±0.31)E+13 




31 


18.2 


(2.10±0.09)E+12 


>6.80E+13 


(1.37±0.09)E+13 


(4.80±0.27)E+12 


(8.53±0.87)E+14 


(3.85±0.16)E+12 


70 


8.7 


(7.76±0.06)E+11 


>6.78E+14 


(3.86±0.64)E+12 


(1.85±0.19)E+12 


(3.72±0.64)E+14 


(8.16±1.09)E+11 


103 


3.1 




>8.11E+13 






(4.13±0.89)E+12 




159 


4.9 


(9.96±3.31)E+10 


>6.64E+13 






(2.25±0.19)E+13 





Table 4. Same as table 3, but for the z=0.91 10 absorber in Q0826-2230. 



Vel 


b eff 


Mgl 


Mg II 


Can 


Fell 


Znll 


-146 


9.4 




(9.77±1.72)E+11 








-126 


11.6 




(4.22±0.28)E+12 




(2.94±0.74)E+12 




-65 


13.1 


(1.08±0.05)E+12 


>1.94E+13 


(3.37±0.19)E+11 


(1.37±0.10)E+13 


(1.82±0.19)E+12 


-37 


14.7 




(2.86±0.26)E+12 








9 


13.8 


(1.42±0.30)E+11 


>1.51E+13 


(8.63±1.83)E+10 


(1.14±0.09)E+13 


(4.26±1.74)E+11 


45 


12.0 




(3.07±0.23)E+12 


(5.38±1.71)E+10 






144 


15.4 




>1.21E+13 


(8.00±1.92)E+10 


(8.98±0.90)E+12 





lines with S/N ~ 20 in the region, so we give the metallicity as an 
upper limit. This system has a low metallicity, with [Zn/H]<-0.98, 
[Fe/H]=-1.28, and [Mn/H]=-1.44 although both Mn and Fe are of- 
ten depleted in the Galaxy. 

System B shows a more complicated velocity structure with 
7 components needed to fit the profile. This system has a larger 
velocity width, spanning ~334 km s . Again, the Al II A 1670 
line was not in the wavelength region, but Al III A A 1854, 1862 
lines were detected. A S/N ~ 25 was reached in the region of Zn 
II A 2026 which was detected at ~ 5a with an equivalent width 
W res t = 41 mA . Profile fitting of the line gives log Nz n n = 
12.36 and [Zn/H]=+0.25. We were able to constrain the ionization 
parameter for this system to be log U=-3.70 based on the observed 
Al III to Fe II ratio (see § 5.3 for a more detailed discussion of the 
method). The Ca II AA 3933, 3969 lines were also detected in this 
system with log Nc a n=12.26. Mn was not detected in this system 
with [Mn/Fe]<-1.24. Tables 5 and 6 show the results of the Voigt 
profile decomposition for systems A and B respectively. 

4.4 Q1010-0047 (z em = 1.671) 

This sys tem has a confirmed sub-DLA with log Nh i=1 9.81 at z = 
1.3270 jNestorll2004l : iRao, Turnshek, fe Nestoj|2006|) . This sys- 
tem was fit with a 8 component model. The profile spanned ~ 308 
km s _1 in total. We detect strong lines of Mg II, Al II, Al III, and Fe 
II in this system but no detection of Zn II A 2026 in the region with 
S/N ~ 30. Table 7 shows the results of the profile fitting for each 
component. From the 3a limiting equivalent width of W rest = 9 
mA , this system has log N Zn n < 11.67 and [Zn/H]<-0.75. We 
detect strong Si II AA 1526, 1808 lines with log N Si n =15.02. This 
system also shows apparent a-enhancement with [Si/Fe]=0.42. Al- 
though depletion may be important, the metallicities based on Si 




Radiol Velocity (km s~') Rodiol Velocity (krr s" 1 ) 



Figure 2. The same as fig. 1 , but for the z=0.91 10 system in the spectrum of 
Q0826-2230 

4.3 Q1009-0026 (z em = 1.244) 

This system has two confirmed sub-DLAs with log Nh i=20.20 at 
z = 0. 8426 (system A) and log N H i=19.48 at z = .8866 (sys- 
tem B) jNestor 112004] : Ir ao. Turnshek. & Nestodl2006h . System A 
shows a relatively simple velocity structure with only 4 compo- 
nents needed to fit the observed profile. The profile for system A 
spans a total of ~94 km s _1 . For this system, the Al II line was 
below the wavelength region covered by the spectrograph, but the 
Al III A A 1854, 1862 lines were observed although the dominant 
species is likely Al II. This system also shows fairly strong Mn II 
AA 2576, 2594, 2606 lines. We did not detect any Zn II AA 2026 
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Table 5. Same as table 3, but for the z=0.8426 absorber in Q1009-0026. 



Vel 




Mgl 


Mg II 


Al III 


Mn II 


Fell 


-28 


8.7 




(1.95±0.21)E+12 








-13 


5.6 




(5.74±0.66)E+12 






(1.53±0.32)E+12 


15 


16.7 


(2.60±0.64)E+11 


>2.83E+13 


(1.16±0.28)E+12 


(4.89±1.11)E+11 




35 


7.2 


(2.73±0.50)E+11 


>1.65e+14 


(4.20±0.33)E+12 


(1.34±0.08)E+12 


(2.47±0.23)E+14 



Table 6. Same as table 3, but for the z=0.8866 absorber in Q1009-0026. 



Vel 


b eff 


Mgl 


Mg II 


Aim 


Call 


Fell 


Znll 



-121 


4.6 




-66 


15.6 


(6 


-13 


23.9 


(1 


70 


5.2 


(1 


87 


3.6 




169 


4.4 




182 


5.1 





.22±0.65)E+11 
.72±0.09)E+12 
.99±0.48)E+11 



(4.47±0.49)E+12 
>7.77E+13 
>2.02E+14 
>2.49E+13 
>1.32E+13 
(1.99±0.31)E+12 
(9.47±1.23)E+12 



(2.51±0.24)E+12 
(9.99±0.36)E+12 
(6.71±2.13)E+11 



(4.58±0.65)E+11 
(1.35±0.08)E+12 



(4.69±0.57)E+12 
(3.22±0.11)E+13 

(3.31±0.31)E+14 (2.28±0.26)E+12 
(2.40±0.52)E+12 
(8.16±0.79)E+12 
(1.93±0.57)E+12 




and Fe are [Si/H]=-0.33 and [Fe/H]=-0.75. As there were detections 
of both Al II and Al III in this system, we were able to constrain 
the ionization parameter U for this system to be log U<-4. 15. See 
§ 5.4 and § 5.5 for a more complete discussion of the Cloudy mod- 
eling and ionization corrections. Figure 5 gives velocity plots for 
some of the lines of interest in this system. The feature at ~-300 
km s _1 in the plot of Al III A 1862 is unidentified. There are several 
other systems seen in the spectrum of this QSO, but no known lines 
match the observed wavelength of this feature. 

4.5 Q1224+0037 (z em = 1.482) 

This QSO has a DLA with log N H i=20.88 at z abs = 1.2346 (sys- 
tem A) and a sub-DLA with log N H i=20.00 at Zq be = 1.2665 
(system B) jNestorbOollRao. Turnshek. fe Nestoi]|2006t) . System 
A shows a relatively simple velocity structure with 5 components 
needed to fit the observed profile. The profile spans ~129 km s - . 
Strong Mg II AA 2796, 2803 and Fe II AA 2344, 2374, 2382 lines 



were detected in this system. Zn II AA 2026 or 2062 lines were not, 
however, detected in this system at S/N ~ 18 in the region, with 
[Zn/H]<-1.62. Strong Al II A 1670 absorption was seen in this sys- 
tem along with weaker Al III A A 1854, 1862 lines. Strong Si II 
AA 1526, 1808 were also seen with log Nsi n=15.10. The weaker 
Fe II AA 2249, 2260 lines were not detected so only a lower limit 
could be placed on log Nf c n > 15.11 based on the saturated Fe II 
AA 2344, 2374, 2382 lines. The metallicities based on Si, Fe, and 
Zn are [Si/H]=-1.32 [Zn/H]<-1.62, and [Fe/H]>-1.24. The Cr II 
AA 2056, 2062 lines were observed in this system with rest frame 
equivalent widths of 47 and 30 mA respectively. The Zn II A 2026 
line was not however detected. Table 8 gives the results of the pro- 
file fitting for each of the components. Figure 6 shows velocity plots 
for several of the lines of interest in this system. 

System B shows a more complicated velocity structure with 
a total of 8 components needed to fit the observed profile. This 
system also shows a larger velocity width of ~253 km s~ . Mg 
II A A 2796, 2803 lines were detected along with strong Fe II A A 
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Table 7. Same as table 3, but for the * = 1.3270 absorber in Q1010-0047. 



Vel 


b eff 


Mgl 


Mgll 


Al II 


Al III 


Sill 


Fell 


-125 


4.9 


(1.53±0.46)E+11 


(8.63±1.59)E+11 










-85 


16.3 


(7.44±0.73)E+1 1 


>3.48E+13 


>7.87E+12 


(4.06±0.37)E+12 


(2.82±0.47)E+14 


(8.72±0.47)E+13 


-55 


4.2 




>4.74E+14 


>5.45E+12 






(6.21±1.61)E+12 


-35 


10.2 




>2.02E+13 


>3.55E+12 


(2.91±0.37)E+12 




(1.65±0.18)E+13 


-14 


3.9 


(3.58±0.62)E+11 


>4.21E+14 


>4.96E+12 


(2.97±0.48)E+12 




(1.56±0.26)E+13 


4 


7.0 




>3.47E+14 


>1.64E+13 


(2.83±0.37)E+12 


(1.66±0.41)E+14 


(7.19±0.73)E+13 


40 


19.9 


(1.62±0.10)E+12 


>1.81E+14 


>3.66E+13 


(6.41±0.43)E+12 


(4.57±0.51)E+14 


(1.31±0.06)E+14 


100 


15.6 




(3.83±0.26)E+12 






(1.43±0.43)E+14 


(7.43±1.43)E+12 


158 


19.2 




(4.18±0.26)E+12 
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Figure 6. The same as fig. 1, but for the z=1.2346 system in the spectrum of 
Figure 5. The same as fig. 1, but for the z=1.3270 system in the spectrum of Q 1224+0037 

QI010+0047 



2344, 2374, 2382 lines. This system also has a saturated Al II A 
1670 line with log N A i n > 13.72 and unsaturated Al III AA 1854, 
1862 lines that were fit simultaneously to give log Nai hi = 12.98. 
The Si II A 1526 line appears to be blended with another feature, 
as can be seen in figure 7. We therefore give an upper limit on 
the column density Nsi ii < 14.30 based on the non-detection of 
the Si II A 1808 line. There also may be some blending near the 
v — —104 km s _1 component in the Fe II A 2600 due to the lack 
of this component in the other Fe II lines. No Zn II AA 2026, 2062 
lines were seen with S/N ~ 20 in the region. Table 9 gives the 
results of the profile fitting for system B. Velocity plots of lines of 
interest are shown in figure 7. 



4.6 Q2331+0038 (z em = 1.714) 

This QSO harbors a sub-DLA w ith log N H i=20.00 dNestor 1120041 : 
iRao. Turnshek. fe Nestor] [2006). This system exhibits a complex 
velocity structure needing 9 components to fit the profile. The to- 
tal profile spans more than 558 km s -1 . The 2 components at 
v — —255 and v = —204 km s _1 appear mainly in the strong 
Mg II AA 2796, 2803 and the stronger Fe II transitions. We de- 
tect Zn II A 2026 in this system at ~ 3a with log Nz n ii = 12.12 
and [Zn/H]=-0.51. The Zn II is found solely in the v — 45 km 
s _1 component, which is the strongest component in Mg I A 2852 
and the Fe II lines. The contribution of Mg I to the blended Zn II 
+ Mg I A 2026 line is small, W ~ 300 mA for the Mg I A 2852 




Rodiol Velocity (km s" 1 ) Rodiol Velocity (krr s" 1 ) 



Figure 7. The same as fig. 1, but for the z=1.2665 system in the spectrum of 
Q 1224+0037 



line, well below the level discussed in § 4. 1 where saturation may 
become important. There appears to be an interloper at v ~ —170 
km s^of the Al III A 1854 line. Table 10 summarizes the results 
of the profile fitting analysis for this system, and figure 8 shows the 
velocity plots for several lines of interest. 
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Table 8. Same as table 3, but for the z=1.2346 DLA in Q1224+0037. 



Vel 


b eff 


Mgl 


Mg II 


AID 


Aim 


Sill 


CrII 


Fell 


-60 


8.2 




(5.85±1.12)E+13 


>2.14E+12 


(7.01±2.02)E+11 






(3.69±1.12)E+13 


-32 


10.8 


(5.14±0.64)E+11 


>9.68E+13 


>2.04E+13 


(4.97±0.34)E+12 


(8.45±0.94)E+14 


(9.27±1.40)E+12 


>3.69E+14 


-9 


8.0 


(8.93±0.84)E+1 1 


>3.73E+14 


>1.50E+13 


(3.29±0.29)E+12 


(3.19±0.75)E+14 


(3.77±1.19)E+12 


>8.51E+14 


16 


10.4 


(7.63±0.74)E+11 


>1.82E+14 


>2.14E+12 


(4.44±0.32)E+12 


(9.16±2.85)E+13 






37 


8.2 




(8.12±1.24)E+12 


>1.69E+12 











Table 9. Same as table 3, but for the z=1.2665 absorber in Q1224+0037. 



Vel 


h eff 


Mgl 


Mgll 


AID 


Al III 


Fell 


-136 


7.8 




>8.95E+13 


>7.11E+12 


(1.00±0.28)E+12 


(4.95±0.94)E+13 


-104 


10.4 




>6.62E+12 


>6.46E+11 






-76 


10.5 




>1.13E+14 


>4.21E+12 


(2.15±0.33)E+12 


(1.07±0.34)E+13 


-26 


5.1 




>3.48E+13 


>1.49E+12 









8.1 


(3.96±0.99)E+11 


>1.23E+14 


>7.16E+12 


(2.75±0.36)E+12 


(5.01±0.96)E+13 


28 


7.9 


(5.93±1.14)E+11 


>1.90E+14 




(7.72±2.71)E+11 


(1.45±0.30)E+14 


57 


7.9 




>8.54E+13 


>1.06E+13 


(1.47±0.30)E+12 


(6.15±1.17)E+13 


92 


8.1 




>3.12E+14 


>2.09E+13 


(1.45±0.30)E+12 


(2.82±0.56)E+13 




Radiol Velocity (km s"') Radiol Velocity (krr s H ) 



Figure 8. The same as fig. 1 , but for the z=l . 1414 system in the spectrum of 
Q233 1+0038 

5 RESULTS 

5.1 Total Column Densities 

In table 1 1 we summarize the total column densities for the differ- 
ent species for all of the absorbers in this sample. Cells with en- 
tries "• • ■" have undetermined column densities because of lack of 
wavelength coverage, being blended with Lyman-Q forest clouds, 
blending with atmospheric absorption absorption bands, or very 
high noise levels due to the inefficiencies of the spectrograph at 
the wavelength extremes. Several systems showed significant sat- 
uration in the lines of Mg II, Al II, and Fe II so the values listed 
are given as lower limits. In general the column densities derived 
from the AOD method are slightly less than the ones derived from 
the profile fitting analysis for unsaturated lines, and differ by up 
to ~ 0.8 dex in the saturated M g II lines. For stronger absorption 
lines, ISavagefeSembacrJdl99lh showed that the instrumental res- 



olution affects the column density derived vi a the apparent opti- 
cal de pth method. For fully resolved profiles, ISavage fe Sembach 
dl99 if) showed that the apparent column density was very close to 
the "true" column density, whereas for unresolved profiles, even 
in moderately strong lines there was significant departure from the 
"true" column density that indicate hidden absorption in the pro- 
file. At the ~20 km s~ resolution obtained with the settings used 
on MIKE, there are likely unresolved components and therefore 
departures from the true column density by those given from the 
apparent optical depth method. 



In table 12 we give the metallicity of these systems from 
profile fitting, and also the abundance ratio [Fe/Zn], which is of- 
ten taken as an indicator of dust depletion. We also give [Si/Fe], 
[Ca/Fe], [Cr/Fe], and [Mn/Fe] for these systems. The ratios of col- 
umn densities of the adjacent ions Mg II / Mg I and Al II / Al III are 
given, as well as Mg II / Al III and Fe II / Al III. For every absorber 
in this sample, there is an under-abundance of Mn relative to Fe. 
This could be a sign of the even-odd effect, although dust depletion 
can confuse this. There also appears to be a trend similar to that of 
DLA systems of higher dep l etion of Fe with increasing metallicity 
dProchaska fe Wolfd 12002!; iLedoux. Petitiean. fe Srianand 12001 
Irlerbert-Fort et al.l 120061 ; iMeiring et al.l |200(j1) seen in these sub- 
DLA systems. Only the z = 1.3270 sub-DLA in Q1010-0047 
shows any signs of a-enhancement with [Si/Fe]=+0.42 although 
dust depletion is also a possibility. It is also worth noting that this 
system shows a low metallicity based on Zn ([Zn/H]<-0.75), but 
higher [Fe/H] and [Si/H], This would suggest that, at least for this 
system, comparing Si and Zn metallicities may give vastly differ- 
ent results. This type of non-standard abundance may not be unex- 
pected, the zones sampled by the background QSOs are small due 
to the small size of the QSO beam, and multiple regions probed 
within a galaxy can show up as different components. Similarly, 
contamination by a small number of supernovae may produce odd 
abundance patterns due to the relatively young, unmixed gas. 
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Table 10. Same as table 3, but for the z=1.1414 absorber in Q2331+0038. 



Vel 


h eff 


Mgl 


Mgll 


A1II 


Aim 


Fell 


Znll 


-255 


10.7 




(7.09±0.76)E+12 










-204 


15.1 




>4.43E+13 


(3.76±0.67)E+12 


(1.61±0.34)E+12 


(1.42±0.22)E+13 




20 


10.4 


(5.77±0.83)E+11 


>2.62E+14 


>8.24E+12 


(1.41±0.31)E+12 


(6.48±0.75)E+13 




45 


16.5 


(1.93±0.13)E+12 


>1.93E+13 


>7.31E+12 


(3.21±0.42)E+12 


(1.02±0.09)E+14 


(1.32±0.32)E+12 


80 


13.8 


(4.50±0.70)E+11 


>2.21E+14 


>4.74E+12 


(1.52±0.33)E+12 


(4.03±0.39)E+13 




186 


11.4 




>3.44E+13 


(1.26±0.40)E+12 


(9.27±2.95)E+11 


(1.16±0.20)E+13 




227 


7.2 




>6.65E+12 










264 


8.3 




>1.35E+13 


(1.13±0.41)E+12 


(1.74±0.31)E+12 


(5.12±0.16)E+12 




299 


5.0 




(2.54±0.49)E+12 











Table 11. Total column densities from the absorbers in this sample. 



QSO 




log N H i 

cm~ 2 


log N Mg i 
cm" 2 


logN Mg II 
cm" 2 


l°g N A1 II 

cm- 2 


log N A1III 

cm- 2 


logN Si n 
cm- 2 


log N Ca u 
cm -2 


logNcrll 
cm -2 


log N M n II 
cm - 2 


log Np e n 
cm -2 


log N Zn ii 
cm -2 


Q0354-2724 


1.4051 


20.18±0.15 


12.72±0.02 


> 15.08 










13.35±0.06 


12.90±0.06 


15.15±0.05 


12.73±0.03 


AOD 






12.70±0.01 


> 14.39 










13.25±0.04 


12.82±0.03 


15.03±0.01 


12.72±0.03 


Q0826-2230 


0.9110 


19.04±0.04 


12.09±0.03 


> 13.76 




< 11.72 


< 14.22 


I1.75±0.06 


< 12.11 


<11.37 


13.57±0.04 


I2.35±0.07 


AOD 






12.06±0.02 


> 13.71 








I1.42±0.04 






13.43±0.12 


I2.41±0.05 


01009-0026 


0.8426 


20.20±0.06 


U.73±0.09 


> 14.30 




12.73±0.05 






< 12.41 


12.26±0.05 


14.39±0.04 


<11.85 


AOD 






1I.80±0.05 


> 13.87 




12.74±0.03 








12.28±0.04 


14.37±0.03 




Q 1009-0026 


0.8866 


19.48±0.05 


12.41±0.04 


> 14.52 




13.12±0.03 


< 14.26 


12.26±0.04 


< 12. 11 


<11.37 


14.58±0.04 


I2.36±0.04 


AOD 






12.43±0.02 


> 14.32 




13.00±0.02 




12.10±0.04 






14.33±0.06 


12.38±0.04 


Q1010-0047 


1.3270 


19.81±0.05 


12.46±0.02 


> 15.17 


> 13.87 


13.28±0.02 


15.02±0.02 




< 12.37 


<11.85 


14.53±0.03 


<I1.69 


AOD 






12.49±0.02 


> 14.26 


> 13.72 


13.29±0.02 


14.86±0.15 








14.50±0.02 




Q 1224+0037 


1.2346 


20.88±0.05 


12.34±0.04 


> 14.86 


> 13.62 


13.13±0.04 


15.10±0.07 




13.12±0.09 


< 12. 14 


>15.U 


<11.89 


AOD 






12.31 ±0.03 


> 14.30 


> 13.59 


13.07±0.O3 


15.01 ±0.05 




13.14±0.05 




> 14.99 




Q 1224+0037 


1.2665 


20.00±0.07 


12.00±0.10 


> 14.97 


> 13.72 


12.98±0.08 


< 14.30 




< 12.37 


<12.16 


14.54±0.09 


< 11.85 


AOD 






12.21±0.08 


> 14.25 


> 13.56 


12.74±0.14 










14.36±0.04 




Q2331+0038 


1.1414 


20.00±0.05 


12.48±0.05 


> 14.79 


> 13.42 


13.02±0.09 


< 14.33 




< 12.37 


<11.79 


14.38±0.03 


12.12±0.11 


AOD 






12.47±0.03 


> 14.39 


> 13.33 


12.86±0.14 










14.44±0.05 


12.22±0.09 



5.2 The Al III to Al II Ratio and Ionized Gas 

Based on the ionization potentials of the Al° and Al + ions (5.99 
and 18.83 eV respectively), one would expect that the majority of 
the Al seen in DLA and sub-DLA systems would be in the first 
ionization stage (Al + ), and that the second ionization stage (Al ++ ) 
would be all but absent due to the shielding of photons with en- 
ergies greater than the ionization potential of hydrogen (13.6 eV). 
The Al III A A 1854, 1862 lines are nonetheless often seen in DLA 
spectra although the majority is always in Al II. As is seen in most 
DLAs, the Al III A A 1854, 1862 lines of the systems studied here 
have very similar component structures to the lower ionization state 
lines such as Fe II, Mg II. Higher ionization state transitions such as 
C IV or S IV often show very different component structures than 
the low ion transitions, implying pockets of local, highly ionized 
gas within the DLA/sub-DLA, or possibly galactic halos. The sim- 
ilar structure of the Al III and low ion profiles suggests that these 
systems are physically connected. 

The Al III / Al II ratios have in the pas t been used to inves- 
tigate the ionization states of QSO absorbers. I Vladilo et al.l l200lh 
observed an anti-correlation between the Al III / Al II ratio and 
Nh i based on a sample of 20 DLAs, where the Al II column den- 
sity was estimated from Si II in most cases. They suggested that the 
higher Nh i systems also had lesser degrees of ionization. It should 
be noted ho wever that there ma y be issues with using Si as a surro- 
gate for Al dBarker et alii 19841) . iDessauges-Zavadskv et ail J2003h 
noted however that when the data were extended into the column 
density region of sub-DLAs, the trend seemed to disappear. In fig- 
ure 9, we show the Al III / Al II ratio for the absorbers from this 



work, Dessauges- Zavad skv et al] J2003). and Vladilo et al. N2001I) . 
As can be seen from the original data of I Vladilo et alj d200ll) (blue 
squares, black diamonds, and purple triangle), the before men- 
tioned trend does seem to hold for these data. It can also be seen in 
figure 9 that the z > 1.8 sub-DLAs of Dessauges-Zavadsk v et al.l 
(2003) and the z < 1.5 absorbers seen here show similar Al III / 
Al II ratios, suggesting little evolution of this parameter although 
the sample sizes are still small. With one exception of 30 cases in 
figure 9, whi ch covers 19<log N H i<21.7, N A i i ii/Nai n <-0.25. 
Models from lPet itiean. Rau ch. fc CarswelH i 19941) also showed lit- 
tle change in the Al III/A 1 II rat io especially when 17 < logNu 1 < 
20. Similarlv. lYork et al] ( l2006t) found the ratio to change very little 
over a wide range of reddenings corresponding to the same range 
in N H i. 

5.3 Cloudy Modeling and Ionization Corrections 

It is generally assumed for DLAs that the gas has a very low ion- 
ization fraction, and that most metal ions are mainly in the first ion- 
ization state. This is because of the self-shielding of photons with 
hu > 13.6 eV due to the large cross section for ionization of H 
at these energies. Therefore, the abundances reported for QSO ab- 
sorbers are typically derived from lines originating from ions and 
may not reflect the true abundances. The ionization correction fac- 
tor, defined here as 

e = [X/H] total - [X+/H ] 

where the total column densities include contributions from all ion- 
ization stages, has been investigated for DLAs by several groups 
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Table 12. Abundances for the absorbers in this sample. 



QSO 




[Zn/H] 


[Fe/H] 


[Fe/Zn] 


[Si/Fe] 


[Ca/Fc] 


[Cr/Fc] 


[Mn/Fe] 


Al III 1 Mil 


Mg II/Mgl 


Mg II / Al III 


Fell/ Aim 


[X/Y] e 




-7.37 


-4.53 


+2.85 


+0.07 


-1.13 


-1.82 


-1.97 










Q0354-2724 


1.4051 


-0.08 ±0.1 6 


-0.50±0.16 


-0.43±0.06 






0.02±0.07 


-0.35 ±0.08 




>2.36 






Q0826-2230 


0.9110 


+0.68±0.08 


-0.94±0.06 


-1.63±0.08 


<0.58 


-0.69±0.07 


<+0.36 


<-0.23 




>1.67 


>2.04 


>1.85 


01009-0026 


0.8426 


<-0.98 


-1.28±0.07 


>-0.31 






<-0.16 


-0.16±0.06 




>2.57 


> 1.57 


1.66±0.06 




0.8866 


+0.25±0.06 


-0.37±0.06 


-0.63±0.06 


<-0.39 


-1.19±0.06 


<-0.65 


<-1.24 




>2.11 


>1.40 


1.46±0.05 


Q1010-0047 


1.3270 


<-0.75 


-0.75±0.06 


>0.01 


0.42±0.07 




<-0.34 


<-0.71 


<-0.71 


>2.71 


> 1.69 


1.25 ±0.05 


Q1224+0037 


1.2346 


<-1.62 


>-1.24 


>0.37 


<-0.08 




<-0.17 




<-0.49 


>2.52 


> 1.73 


>1.98 




1.2665 


<-0.78 


-0.93 ±0.11 


>-0.16 


<-0.31 




<-0.35 


<-0.41 


<-0.74 


>2.97 


> 1.99 


1.56±0.12 


Q2331+0038 


1.1414 


-0.51±0.12 


-1.09±0.06 


-0.59±0.11 


<-0.12 




<-0.19 


<-0.62 


<-0.40 


>2.31 


> 1.77 


1.36±0.09 



This work Upper Limits (z<l .5) 
D-Z {2003) Detections ( 1.8<z<4.3) 
D-Z (2003) Upper Limits (1.8<z<4.3) - 
Vladilo (2001) Estimations (1.6<z<4.5) 
Vladilo (2001) Upper Limit (z=2.8) 
Vladilo (2001) Al ++ / AI + (z-2.2) 



ih^ i 



-L -r- ib 



20.5 
log N(H I) 



Figure 9. The Al III / Al II ratio for th e abso rbers in this sample 
(striped triangles), [P essauges-Zavad skv et al"] j2003l) (solid triangles), and 
IVladilo et aifbOOlT ( blue s quares, purple triangle, and black diamonds). 
For the IVladilo et all 120011) systems, the blue squares represent systems 
where Al II was estimated from Si II, the black diamond points are for 
systems with true Al III / Al II values, and the purple triangle is an upper 
limit based on Al III / Al II. All of the absorbers from this work with Al II 
measurements showed significant saturation so only lower limits could be 
placed on their column densities. 



( Howk <fe Sembacbi 1 19991 ; IVladilo et all 1200 ll : iProchaska et~ai1 
l2002h . Their conclusions were all roughly the same; that the ion- 
ization corrections for most elements were <0.2 dex in most cases. 
Due to the smaller amount of HI in the sub-DLA systems, one 
mig ht expect that they would show a g reater amount of ionization. 

iDessauges-Zavadskv et alj J2003f) showed the ionization cor- 
rections for their sample of sub-DLA at 2 > 1.8 were, in gen- 
eral small, with <0.2 dex corrections based on predictions of 
the Cloudy software package. The only absorber in their sam- 
ple that showed significant signs of ionization effects was the 
log N H i=19.37 sub-DLA in the spectrum of Q2155+1358 at 
Zaba = 3.565. The absorbers in the spectra of Q0826-2230 with 
log N H 1 = 19.04, Q1009-0026 with log N H i=19.48, and Q1010- 
0047 with log Nh i = 19.81 have the lowest column densities in 
our sample, and thus might be expected to require some ionization 
corrections. 

To investigate further, we used the Cloudy software package to 
compute photo-ionization models assuming local thermodynamic 
equilibrium (LTE). Calculations were performed with version 



C06.02.b of Cloudy, last described by iFerland etai] dl998h . Grids 
of Cloudy models were computed assuming that the spectrum of 
ionizing radiation striking the cloud followed the form of the extra- 
galac t ic UV background of Haardt & Madau dHaardt fcl vladau 
1 19961 : iMadau. Haardt. fc Reesll 19991) at the appropriate redshift of 
the absorber, plus the model stellar atmosphere of Kurucz with a 
temperature of 30,000K to simulate a radiation field produced via 
an O/B-type star. Plots of both types of spectra can be found in 
Hazy, the documentation for Cloudy. The extra-galactic spectrum 
of Haardt & Madau is a harder spectrum than the Kurucz spectrum, 
which produces few photons with hi/ >2 Ryd. In fact, the Haardt 
& Madau spectrum produces a total of ~ 21.13 times as many 
photons with hv > 13.6 eV than the Kurucz spectrum. From 1-2 
Ryd though, the total ionizing photons is almost identical for both 
spectra. It has recently been shown that the contribution from local 
sources to the ionization of DLA systems is likely non-neg ligible 
compared to the ionizing background radiation field ( Schave 2006). 
We include the radiation field of the Kurucz stellar spectrum in 
these Cloudy simulations, but do not allow that intensity to vary 
relative to QSOs. In addition, our Cloudy models have included 
both cosmic ray and cosmic microwave backgrounds. These simu- 
lations however do not include radiation from local shocks caused 
by supernovae, or compact sources such as white dwarfs or com- 
pact binary systems, all of which may contribute significantly to the 
ionizing radiation field. Explorations into contributions of various 
galactic sources is beyond the scope of this work. 

For each of the grids of models, the ionization parameter de- 
fined by 



U 



<t> 9 



nii 



(where $912 is the flux of radiation with hv > 13.6 eV) was in- 
creased from log U=-6.0 to 0. Each model was tailored to match 
the Nh i, metallicity based on Zn II, and redshift of the observed 
system. Due to the saturation of the strong Al II A 1670 transition, 
for each of the three systems where we have detections of Al III AA 
1854, 1862 lines only an upper limit can be placed on the ioniza- 
tion parameter U. Also, we are not at all able to constrain the shape 
of the ionizing spectrum that is illuminating the cloud. The pos- 
sibilities presented here are but a myriad of other possible shapes 
of the spectrum of ionizing radiation. Nonetheless, using standard 
assumptions, outlined above, some conclusions can be drawn. 

Cloudy models were calculated for all of the systems ob- 
served. We present the results of the most interesting systems here. 
With log N H i = 19.04 the weak sub-DLA toward Q0826-2230 may 
be expected to require ionization corrections. The only element ob- 
served in this system with multiple ionization stages was Mg. Due 
to the saturation of the Mg II AA 2796, 2803 lines only a lower 
limit could be placed on the Mg + to Mg° ratio, and thus a lower 
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limit on the ionization parameter U. However, the Mg II A 2803 
line was not optically thick even in its strongest components, so the 
column density given in table 12 may not be far from true. To see if 
the true column density could be significantly higher, we increased 
the column densities in the components by a factor of 10, at which 
point there was significant deviation between the observed and the- 
oretical profiles. Also, the column densities derived from the AOD 
method and profile fitting agreed to within 0.05 dex. Nonetheless, 
as can be seen from figure 10, the ionization corrections to Zn + are 
<0.2 dex even for log U<-3.0 which is the typical upper limit for 
other sub-DLA systems observed. We have therefore not included 
any ionization correction factor into the quoted Zn abundance in 
table 13. 

In figure 10, the horizontal lines in the center panels indicates 
the observed limits or ratios of the Mg II to Mg I, Al III to Al II, and 
Al III to Fe II column densities. The vertical lines are therefore the 
corresponding ionization parameters based on these values. Based 
on the Mg II to Mg I ratio in the z aha = 0.9110 sub-DLA in the 
spectrum of Q0826-2230, the ionization parameter is a lower limit 
due to the saturation of the Mg II lines. 

The log Nh i = 19.48 system in Q1009-0026 may also require 
ionization corrections. The Al II A 1670 line was below the wave- 
lengths accessible with the MIKE spectrograph, but we the Al III 
A A 1854, 1862 lines were observed. We were therefore able to con- 
strain the ionization parameter based on the Al III to Fe II ratio, 
which was log(Al III / Fe H)=-1.66, and consequently based upon 
the Cloudy models log U=-3.70. The ionization correction for Zn + 
based on this ionization parameter is small (~ 0.15 dex), so no cor- 
rection corrections were introduced into table 13. 

The log N H i = 19.81 system in the spectrum of Q1010-0047 
is another system with relatively low Nh i that may require ioniza- 
tion corrections. For this system, both the Al II A 1670 and the Al 
III A A 1854, 1862 lines were observed. The Al II A 1670 line is 
saturated, so only an upper limit could be placed on the A1++ to 
Al + ratio. From the Cloudy models, the ionization parameter log 
U<-4. 15. Again, the ionization corrections for Zn + were small for 
values within the range of possible values for U, so no ionization 
correction was added to the values in table 13. 

In general, from results of the Cloudy modeling the ionization 
parameter U is typically small, log U<-3.0 for the systems that had 
detections of the Al II A 1670 and Al III A A 1854, 1862 lines. We 
note, again, that the ionization parameter does depend strongly on 
the shape of the ionizing spectrum. 

It has been mentioned before that the Al III to Al II ratio 
may not be an accurate predictor of the ion ization state of the ga s 
observed in QSO absorption line systems jProchaska et al.ll2002l) . 
One of the reasons for this is the uncertainties in the atomic data for 
Al. It has been suggested that the Al III to Al II ratio may provide 
an estimate of the ionization of the gas in DLAs to first order, but 
better diagnostic indicators of the ionization state of the gas are the 
Fe ++ to Fe + and N + to N° ratios. For the absorbers presented here 
the lines of these ions are well in the UV therefore could not be 
observed with our ground based observations. Therefore the Al III 
to Al II ratio is the only opportunity available to study the ioniza- 
tion state of the gas in these sub-DLA systems. Even if we cannot 
derive the ionization parameter absolutely, it appears that the ion- 
ization of the components does not change over a wide range of 
integrated column densities and redshifts. It should be noted that 
for the systems in which Al II was not observed there is no indi- 
cation that Al ++ is the dominant ion. As can be seen in table 14, 
the Mg II to Al III ratio is similar for systems with and without 
Al II measurements. If Al III were relatively higher in one sys- 
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HM +Kurucz spectrum 
log N H =19.04, z ibi =0.91 10 




logU 



Figure 10. Results of the Cloudy simulations for the z = 0.9110 sub- 
DLA with log N H i=19.04 in the spectrum of Q0826-2230. The upper panel 
shows the column densities for several species of interest. The middle panel 
shows the logarithmic ratio of Mg+ to Mg°, A1++/A1+, and Al++/Fe+. 
Here, the horizontal solid red line indicates the lower limit to the observed 
ratio of the column densities derived from the Mg I A 2852 line, and the Mg 
II AA 2796, 2803 doublet. The vertical solid blue line is then the lower limit 
of the ionization parameter U, based upon the Cloudy modeling. The dashed 
lines show the observed Al III to Fe II ratio (purple), and the corresponding 
upper limit to the ionization parameter (yellow). The lower panel shows the 
ionization correction factor for Zn + in dex. 



tern than the others, the limit of Mg II / Al III would be lower 
than the others, r The use of the Al III / Fe II ratio in determining 
the ionization parameter may introduce some uncertainties. While 
Cloudy assumes solar abundances, differential depletion between 
Al and Fe may lead to an underestimation of the Al III / Fe ratio. 
In the Galaxy, Al i s sometimes seen to be more depleted than Fe 
jBarkeretal1ll984l) . Another source of uncertainty are nucleosyn- 
thetic mismatches between these two elements. The use of column 
densities from adjacent ions such as Al III and Al II, or Mg II / Mg 
I assumes that the ions coexist in the gas from which these lines 
arise. This is almost universally true in the case of Al III and Al II, 
the line profiles are very similar. The line profiles for Mg II and Mg 
I are not always as similar, especially in the components at higher 
radial velocities where Mg II lines are often seen without any com- 
ponents seen in the Mg I profile. The "core" components however 
typically do match very well, as can be seen in figures 1-8. 
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5.4 Sub-DLAs and the velocity dispersion-metallicity 
relationship 

Based on a sample of ~ 53,000 star-forming galaxi es at z~0. 1 ob- 
served in imaging and spectroscopy in the SDSS, Tremonti et al.l 
d2004l) discovered a mass-metallicity relationship for these galax- 
ies. Specifically, they found a correlation between stellar mass and 
metallicity that spans over 3 orders of m agnitude in s t ellar m ass 
and one order of magnitude in metallicity. iNestor et alJ J2003h no- 
ticed a relation between the width of Mg II and the metallicity, 
and suggested that the Mg II line width might be an indicator that 
the gas was in a deep potential well. Evidence has recently been 
provided for the possible existence of a mass metallicity relation- 
ship for DLA absorbers, assuming the velocity width of op tically 
thin lines to be proportional to the mass dLedoux et alj|200q) . As a 
proxy for the stellar mass of these systems, which has been difficult 
to detect, the velocity width is used as an indicator of galaxy mass, 
as it potentially probe s the depth of t he gra vitational potential well 
of the DLA systems. iBouche et all (2006) however find an anti- 
correlation between the Mg II equivalent width, and the estimated 

halo mass based upon an indirect mass indicator . 

Following the analysis of iLedoux et alJ |2006), we per- 
formed an analysis of the apparent optical depth (see e.g., 
ISavage fc Sembacr] (fl996)) of these sub-DLA systems. The appar- 
ent optical depth is defined as 

T(\) a = ln[I (\)/I obs (\)], 

where Io(A) is the continuum level, and I 6 S (A) is the observed 
intensity. Specifically, we defined the beginning of the absorption 
line systems profile as the point where the apparent optical depth of 
the line reached a value of 3a higher than the average noise value in 
the continuum. Similarly, the upper edge was defined as the point 
where the apparent optical depth dropped below the 3 a level of the 
continuum. In figures 1 1 and 12, we give the apparent optical depth 
of the Fe II A 2374 lines vs. the radial velocity widths defined here 
for the systems in this sample. 

In figur e 13 we plot the metal licity vs velocity for the DLA ab- 
sorbers fromlLedoux et al.l J2006t) (open circles); su b-DLAs from 
ILedoux et all f2 006) (o pen squ ares ); sub-DLAs fromlperoux et al.l 
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Figure 11. Apparent Optical Depth plots for the Fe II A 2374 line for 
(clockwise from top left); the z abs = 1.4051 system in Q0354-2724, the 
z abs = 0.9110 system in Q0826-2230, the z abs = 0.8866 system in 
Q1009-0026, and the z aba = 0.8426 system in Q1009-0026. The total 
velocity width for the system, defined in § 5.3, is given above the plot. 



components at higher radial velocities that are often not seen in the 
intrinsically weaker features. 

As can be seen from figure 15, the systems from this sample 
are in agreement with the general trend of the data. More interest- 
ing though is the fact that the sub-DLAs from this sample and from 
Prochaska et al. (2006); Peroux et a i] j2006ah are all at the highest 
end of the metallicity range. Indeed, the only points that lie above 
solar metallicity are those from sub-DLAs. Is this possible evidence 
that sub-DLAs are more representative of massive sp iral or ellip- 
tical g alaxies than their DLA counterparts? See also iKhare et al.l 
(2006) for a discussion of the nature of QSO absorbers. 

Because the absorbers chosen in our sample were chosen 
partly because of their strong Mg II or Fe II lines, it is possible 
that these systems would be disposed to lie in the upper right cor- 
ner of figure 13, with high metallicities and larger velocity widths. 
Ledoux et all II2UU6I) (open squ ares ); sub-DLAs iromlFeroux et all It was, however, show n that Wm s ii does not correlate with [Zn/H] 
1 2006al) . IPeroux et al.l (I2006bh. fc IProchaska et alJ J2006I) (filled jKulkarni et"ai]|2007l) based on individual measurements. The iso- 
squares); DLAs from lPeroux e t al. ( 2006al) fc lPeroux et aTT l2"b06b) lation of points in this region is therefore a real effect, 
(filled circles); the Zn detections from this work (filled stars); and 
the up per limits from th is work (downward arrows). In their analy- 
sis, |Eedoux~etal] {200(J) used weaker features, where at most 40% 
of the continuum level was absorbed by the line. Due to the un- 
availability of the spectra for the sub-DLAs in the literature, we 
were unable to measure the velocity width directly, and the value 
plotted is the velocity width given in the ref erence. In the case of 
the two sub-DLAs in Prochask a et al. I J20061) . the velocity width 
was estimated from the plots of the Si II A 1808 and Fe II A 16 08 
lines. For the systems in our sample and IPeroux et al.l j2006al fbh. 
in addition to the Fe II A 2374 line we also measured the velocity 
width of the Mg II A 2796 or A 2803 line, which is also plotted 
for each system in figure 15, with a line connecting the two values. 
We have chosen these two lines for this investigation for the several 
reasons. Firstly, with weaker features, the velocity width becomes 
sensitive to the S/N. Secondly, the Fe II A 2374 line looks much 
like Ca II in the SMC, so ground based work with Ca II in lower 
z systems can be compared directly with the Fe II A 2374 line in 
higher z systems. Lastly, the stronger Mg II A 2796, 2803 lines rep- 
resent a maximum velocity width to the line, and also include the 



6 CONCLUSIONS AND FUTURE WORK 

In this paper, we have presented high-resolution abundance mea- 
surements of 7 sub-DLAs and 1 DLA. Several of these systems 
show moderately high levels of metal enrichment compared to 
DLAs, with three having near solar or super-solar abundances. 
To date, the higher HI column density DLA systems with log 
Nh i > 20.3 have been studied preferentially due to their high 
gas content and therefore possible high metal content. However, 
these systems show low metallicities at any redshift with [Zn/H]~- 
1.0 dex. Of all the (mainly neutral) QSO absorbers studied to date, 
all eight that h ave been observed with super-solar m e tallicities are 
all sub-DLAs jPettini et~aT]|200(j|; IKhare et al.ll2004l; IPeroux et all 
2006a; IProchaska et al.l2006l) . From the results of the Cloudy mod- 
eling presented here, the ionization corrections for these absorbers 
in general app ear to be low, consistent w i th the ionization correc- 
tions found by|5 essauges-Zavadskv et a 0( 120031) for the sub-DLAs 
at higher redshift. We note however that the shape of the ioniz- 
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Figure 12. Same as figure 13, but for the z a i, s = 1.3270 system in Q1010- 
0047, the « ai)s =1.2346 system in Q1224+0037, the z a6s =1.1414 system in 
Q2331+0038, and the z ai)a =1.2665 system in Q1224+0037 
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Figure 13. T he metallicity vs v elocity width for the DLA ab- 
sorbers from iLedoux et al I 120061) (open circl es); sub-DLAs from 
Ledouxe^alj i 2006J) (ope n squares); sub-DLAs fr om lPeroux et all i2006al) . 
Peroux et alf . 2006bl). &. IProchaska et alj 120061) (filled squares); DLAs 
from lPeroux et alj j20Q63) fc lPeroux et all l2006bl) (filled circles); the Zn 
detections from this work (filled stars); and the upper limits from this work 
(d ownward arrows). Fo r t he absorbers in this sa mple and for the absorbers 
in lPeroux et all l2006ah fc lPeroux etal] fe006bl) . we measured the velocity 
width from both the Fe II A 2374 line as seen in figures 13 and 14, and 
also the stronger Mg II AA 2796,2803 lines. These two values are plotted in 
figure 15 with a line connecting the two points. 



rors for all our values are based on formal errors. We cannot quan- 
titatively evaluate possible sources of error from external sources 
outside the range of the assumptions used in the analysis of Nh i 
and the ionization corrections. The matter of Nh i determinations 
is discussed in the appendix, and ionization corrections are dis- 
cussed in § 5.3. The former requires higher resolution, higher S/N 
UV spectra to improve the confidence in the errors, and the latter 
requires a much more detailed knowledge of the radiation field in 
the absorbers as well as the geometry. Nh i determinations could 
potentially be improved when HST is refurbished with UV spec- 
troscopic capability. The ionization corrections are more difficult 
to determine, and will need to be continually improved when more 
is learned about the environment of the absorbers. For the moment, 
we urge caution when interpreting the highest abundances, until 
more cases are verified or these two matters are more fully dealt 
with. 

The errors reported in table 12 ar e based upon the errors 
given bv lRao. Turnshek. fe Nested d2006l) . with no ionization cor- 
rections. Taking the errors for Nhi, conservatively estimated in 
the appendix, our estimates on the ionization correction from the 
cloudy modelling, and the error on Nz n ii> the metallicities for 
the z abs =0.9110 and z abs =0.8866 sub-DLAs in Q0826-2230 and 
Q1009-0026 are +0.68±0. 17 and +0.25±0. 16 respectively with the 
errors added in quadrature. We have however presented arguments 
that the observed trends in ion ratios over large ranges of Nh i make 
our conclusions valid, for the purposes of this paper. 

Although progress has been made recently in the understand- 
ing of sub-DLA absorbers, more observations are needed to fully 
understand the role of sub-DLAs in the metal enri chment his- 
tory of the universe. Aside from the two absorbers in lPeroux et al.l 
our measurements constitute the only available high res- 
ol ution data available for sub-DLAs at z < 1.5. As was pointed out 
in lProchaska et alj d2006t) . these very metal rich sub-DLAs may be 
rare, but even if they are < 1 % of the total population they will con- 
tribute significantly to the mean metallicity because of the much 
larger number density of sub-DLA systems. Clearly, more observa- 
tions are needed to determine how common metal-rich sub-DLAs 
are, to see if sub-DLAs do harbor a large fraction of the missing 
metals, and to extend the mass-metallicity relationship for QSO ab- 
sorbers to include more sub-DLAs. Future HST observations would 
allow us to better determine the ionization corrections using the Fe 
III A 1 122, N II A 1084, and N I AA 1 134, 1 199 lines respectively. 
The two transitions of S III (AA 1012, 1 190) provide a more likely 
unblended source of lines that can be used to check for ionization 
corrections, though S III probes somewhat harder radiation fields 
than Fe III and N II. 



ing spectrum did strongly effect the ionization corrections from the 
Cloudy models, so the true ionization corrections could possibly 
be much larger. We have based our calculations on the radiation 
fields currently thought to be most relevant. Conditions will vary 
from component to component, and from system to system based 
on their respective temperatures and densities, and on the proxim- 
ity or not of X-ray sources, UV stellar sources, and UV emission 
from shocks. The input spectrum is therefore certainly a source of 
uncertainty, but other input parameters ranging from atomic data to 
the morphology assumed by Cloudy in modeling DLA systems are 
likely sources of uncertainties. 

On the matter of the high (i.e. super-solar) abundances, the er- 
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APPENDIX - N H i DETERMINATIONS 

The systems studied in this work all have known Nh i from HST 
spectra. Here we provide plots of the Voigt profiles of the Lyman- 
q transition using the best fit v alues of the column density from 
iRao. Turnshek. fc NestoJ |2006). Due to the low resolution and S/N 
of these UV spectra, only one component was used in the fits. Sev- 
eral factors such as blending with Lyman-a forest lines, contin- 
uum fitting and component placement add difficulties to the pro- 
file fitting with these spectra. These effects are largely unknown 
in each case and formal errors that account for a wide range of 
possibilities cannot be made until there are higher S/N UV spec- 
tra from the refurbished Hubble Space Telescope. We show in 
figures 14 and 15 the Voigt profiles correspondin g to the column 
densities given bv lRao. Turnshek. fc NestoJ fc006l) and convolved 
with the instrumental spread function, superimposed on the archival 
data from HST cycle 6 program 6577, cycle 9 program 8569, 
and cycle 11 program 9382. We note that the normalization, i.e., 
the continuum fit that we defin e may differ from that adopted by 
IRao. Turnshek. fc NestoJ ^20061) . For our continuum fits, a polyno- 
mial typically of order <6 was used, and the absorption line itself 
was excluded from the fitting region. Also over-plotted are pro- 
files with H I column densities smaller and larger by 0.15 dex. For 
the sub-DLAs in our sample, the plots are not in print, and we in- 
clude the DLA for completeness . In some cases, the errors cited by 
IRao. Turnshek. fc NestoJ d2006t) may not cover the revised values 
from higher resolution, higher S/N data, and the plotted range gives 
an idea of how much things may change. However, we find the fits 
and errors in column density from lRao. Turnshek. fc NestoJ ( 2006) 
to be completely satisfactory, and have used them in the body of 
the text. In cases where the profiles with error profiles that seem 
to fit slightly better, there is no systematic sense of the differences: 
in some cases the H I column density is smaller, and in others it is 
larger. This may imply that our values of [Zn/H] may need revision 
later, when higher quality data become available. In particular, the 
system in Q0826-2230 may have a slightly higher value of Nh i 
and therefore slightly lower [Zn/H]. However, the value would still 
be ~2x solar. The conclusions of this paper would not be changed 
by any likely changes in the values of Nh i. 
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